We plot a variable of interest from CMAP in a range of latitude and longitude at a given point of time.
Specify following parameters to retrieve required data from CMAP.
plot_regMap function download the data, and output a list which contain a) regional map plot object: plot_ly/ggplot with “x axis” as “longitude” and “y axis” as “lattitude”; b) corresponding data table.
con <- dbConnect(odbc::odbc(), DSN="CMAP-MSSQL",UID="ArmLab",PWD="ArmLab2018")
# Inpit variable:
table.name = 'tblsst_AVHRR_OI_NRT'
sel.var = 'sst'
range.var <- list()
range.var$lat <- c(10,70)
range.var$lon <- c(-180,-80)
range.var$time <- c('2016-04-30', '2016-04-30')
## Data retrieval
tbl.subset <- getTableData(con, table.name, sel.var, range.var, order.var= c('lat','lon'))
head(tbl.subset)
## # A tibble: 6 x 4
## lat lon time sst
## <dbl> <dbl> <date> <dbl>
## 1 10.1 -80.1 2016-04-30 29.6
## 2 10.1 -80.4 2016-04-30 29.8
## 3 10.1 -80.6 2016-04-30 29.9
## 4 10.1 -80.9 2016-04-30 29.9
## 5 10.1 -81.1 2016-04-30 29.8
## 6 10.1 -81.4 2016-04-30 29.8
## Plot - regional map
# out <- plot_regMap(con, table.name,sel.var,range.var, type = 'ggplot')
out <- plot_regMap(con, table.name,sel.var,range.var, type = 'plotly')
out$plot
dbDisconnect(con)
con <- dbConnect(odbc::odbc(), DSN="CMAP-MSSQL",UID="ArmLab",PWD="ArmLab2018")
table.name <- "tblPisces_NRT"
sel.var <- "Fe"
range.var <- list()
range.var$lat <- c(10,70)
range.var$lon <- c(-180,-80)
range.var$depth <- c(0,0.5)
range.var$time <- c('2016-04-30', '2016-04-30')
## Subset selection:
tbl.subset <- getTableData(con, table.name, sel.var, range.var, order.var= c('lat','lon'))
head(tbl.subset)
## # A tibble: 6 x 5
## lat lon depth time Fe
## <dbl> <dbl> <dbl> <date> <dbl>
## 1 10 -80 0.494 2016-04-30 0.000725
## 2 10 -80.5 0.494 2016-04-30 0.000705
## 3 10 -81 0.494 2016-04-30 0.000658
## 4 10 -81.5 0.494 2016-04-30 0.000635
## 5 10 -82 0.494 2016-04-30 0.000640
## 6 10 -82.5 0.494 2016-04-30 0.000603
## Plot - regional map
# out <- plot_regMap(con, table.name,sel.var,range.var, type = 'ggplot')
out <- plot_regMap(con, table.name,sel.var,range.var, type = 'plotly')
out$plot
dbDisconnect(con)